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(<~^ . Abstract: We obtain ratchet effect in inertial structureless systems in symmetric 

o ; 

\ periodic potentials where the asymmetry comes from the nonuniform friction offered 



by the medium and driven by symmetric periodic forces. In the adiabatic limit the 
calculations are done by extending the matrix continued fraction method and also 
by numerically solving the appropriate Langevin equation. For finite frequency field 
' drive the ratchet effect is obtained only numerically. In the transient time scales the 



O 



X 



system shows dispersionless behaviour as reported earlier when a constant force is 
applied. In the periodic drive case the dispersion behaviour is more complex. In this 



\ brief communication we report some of the results of our work. 
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1^ ■ I. INTRODUCTION 

o 

O \ The phenomenon of obtaining net unidirectional current in a periodic potential with- 

^ ■ out the application of any time averaged external field is termed as ratchet effect 
00 ! 

Q , This is necessarily a non-equilibrium phenomenon and has been investigated extensively 



in systems where damping to periodic motion is large. Also, in majority of the investi- 



gations the potential is considered asymmetric (and hence the name ratchet). In order 

B • n 

to obtain ratchet effect in asymmteric potentials the system is rocked [2j periodically 



(rocked ratchet) or the potential amphtude is changed[3|] dichotomously, periodically or 
randomly, (flashing or fluctuating-potential ratchets). The effect can also be obtained 
if the system is driven periodically but time asymmetrically in such a way that the 
total applied force per period is zerojj]. In such systems the condition for asymmetric 
periodic potential can be waived. All these popular models consider the systems to 
be homogeneous where the friction coefficient is taken to be constant in space and the 
temperature is maintained uniform. In the present work we consider an underdamped 
system and the periodic potential is taken to be symmetric. However, the friction 
coefficient is considered to vary periodically in space similar to the potential but with 
a phase difference. Also, the system is driven periodically symmetrically about zero. 
The overdamped case of the problem has been studied earlier {s]. 



riction coefficient 
dl. Nevertheless, 



Since the system asymmetry is solely due to non-uniformity of 
the expected net particle current in such a system may be weak 
the physical explanation for asymmetric particle current can be easily understood. 
In fig.l the periodic potential V{x) = —sin{x) and the friction coefficient j{x) = 
7o[l — Xsin{x + (p] are plotted together with A = 0.9, phase difference = 0.35 and 
7o = 1. In the figure a period of the potential is shown divided into two halves one with 
friction coefficient small (< 1.0) and the other with large (> 1.0) friction coefficient. 
The effect of (p on this division is clearly visible: The friction is more on the left of the 
potential peak position than on the right. 

The physical explanation for the possibility of obtaining ratchet effect in this in- 
homogeneous system can be given in two ways. Firstly, since in the static situation 
the position probability distribution is independent of the profile of the friction coef- 
ficient, the distribution will be Gaussian-like and symmetric about the minima of the 
potential. However, the situation changes in the dynamic case. Because of the higher 
friction on the left of the potential peak position the particle spends more time there 
during its motion. And hence the effective position probability distribution is skewed 
with higher probability on the left than on the right of the potential peaks as though 
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FIG. 1: Plot of potential V{x) (solid line) and the friction coefficient 7(x) (dashed line). The 
portion of a period of V{x) where 7(x) < 1 is shown hatched and the portion where 7(x) > 1 
is shown by a thick line. 
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a static constant force has been applied |5i] in the left (negative) direction. Thus, one 
would expect a negative average particle current even without the application of net 
external field. Secondly, since the particles tend to spend more time on the left of the 
peak it will absorb more thermal energy from the bath and hence effectively it will 
be at a higher temperature on the left side of the potential peaks. Higher effective 
temperature on the left leads to higher probability of thermally assisted passage from 
the left of the potential peak to the right giving rise to a net current in the right (or 
positive) direction. These two opposing effects in the nonequilibrium situation may 
lead to competition in deciding the direction of the ratchet current depending on the 
temperature. These two competing effects could even lead to ratchet current direction 
reversals as will be revealed by our numerical results. 

The underdamped systems are difficult to study (analytically as well as numerically) 
but provide interesting results. If the particle is initially put at the bottom of one of 
the wells of the potential (not with great initial velocity) it will get out of the well 
due to thermal effect (Kramers-like behaviour). However, once out of the well it will 



roll down from the top to the bottom of the next well of the tilted periodic potential 
without much frictional loss of energy and hence may continue to move along for a long 
distance (more than one period) before it again gets trapped in a distant well. From 
there it will again get out due to thermal effects and the process will keep repeating. 
The distribution of travel distances between two consecutive haltings is itself quite 
interesting to study. When taken an ensemble of such particles they may collectively 
show coherence properties in the intermediate time scales, where position dispersions 
< (Ax(t))^ >=< x{t)^ > — < x{t) >^ become constant in time for a brief period. This 
is how the long time steady state behaviour < {Ax{t))^ >~ t is bridged with the very 
short time behaviour < (Ax(t))^ >~ t'^. This result has been reported recently for a 
constant applied forceQ]. In the periodically driven case, however, the system shows 
more complex behaviour to be described as we proceed. 

Obtaining analytical solutions of the Langevin equation for underdamped inhomo- 
geneous (space-dependent friction) systems subjected to constant uniform force field 
in the steady state is in principle a straightforward extension of the matrix continued 
fraction method (MCFM) initially developed by RiskenQ]. Our extended formalism, 
however, yields stable solutions only in a limited range of the applied force, especially 
for small friction coefficients. It turns out that the range of the applied force where 
MCFM yields sensible results is the one relevant to our discussion. In section II, we 
simply describe how to calculate the mean particle velocity in the steady state. How- 
ever, it has not yet been possible to write down a similar expression for the solution of 
the Langevin equation when the system is subjected to a time varying external field. 
We, therefore, resort to numerical methods. The efficacy of the numerical method is es- 
tablished by comparing the results so obtained with the MCFM result in the adiabatic 
case (Fig. 2). The numerical results along with the MCFM (adiabatic case) results are 
presented in section III. In section IV we discuss the salient features of our work to 
conclude. 



II. THE METHOD 



Consider a particle of mass m moving in a periodic potential V{x) — —Vosin{kx) 
in a medium where the particle experiences friction with coefficient 7(x) = 70 (1 — 
\sin{kx + (p)) with A < 1. The motion can be described by the following Langevin 
equation: 

t 

dt^ dt dx 

where T is the temperature in energy units in terms of fcs, the Boltzmann constant and 

F{t) is the externally apphed field, ^{t) is the fluctuating force term with < ^{t) >= 0, 

and < ^(t)^(t') >= 25{t — t'). In dimensionless units with m = 1, Vq = 1, = 1, the 

Langevin equation reduces (reduced variables denoted by the same symbols) to 

d^x dx I 

^ " "'^^^^ + ''''''' + ^^^^ + V7(^)^eW, (2) 
where 7(0;) = 70(1 — Xsin{x + (p)) and with similar noise statistics for the reduced ^{t) 
as earlier. Our purpose is to calculate x{t) and v{t) by solving the Langevin equation. 
As mentioned earlier, so far no analytical solution could be found even for the average 
velocity of the particle when the external force F{t) is time varying. However, when 
F{t) — Fq, a constant, it is straightforward to extend the MCFM to obtain the drift 
velocity in the steady state situation. The method involves solving the Fokker-Planck 
equation 

=£,,iy (3) 

corresponding to the Langevin equation (2), where the FP operator 

d d d 

C-TV -^7^ + 7o(l - sin{x + 0))Tr^ - {cosx + -Fo)— + 
ox ov ov 

-foT{l-Xsin{x + 4>)) — . (4) 

In the MCFM, the distribution W{x,v,t) is expanded in terms of the function ipn, 

2 °° 

W{x,v,t) = {27rT)^e^J2Cn{x,t)i;n, (5) 

n=0 



where i/jn = ^^/^° ^'i'^ the eigen functions of the operator 'y{x)b^b and 

11)0 = {2n)^T^e^. (6) 

b = VT-^ + b'^ = ~VT-^ + steady state, the Fokker-Planck equation 

reduces to 

n-f{x)Cn{x) + Vn + lDCn+i{x) + V^DCn-i{x) = 0, (7) 

where D = VT-^ and D = \/T^ + . This is a series of recurring equations for 

n = 0,1,2, N. This series is truncated at = which depends on the value of 
7, F, etc and can be quite large. Since the potential as well as the friction coefficient 
are periodic, the mean particle velocity < > is given by dv J^"" dxW{x, v)v which 
turns out to be a constant from the first, n = 0, of the series of the equations (7). 
Solving these simultaneous equations to obtain Cn{x) is not straightforward, however, 
and one may not obtain any stable solution for all values of 70, and F. 

For the time dependent external force fields, one necessarily have to use numeri- 
cal methods to solve the Langevin equation. For homogeneous inertial systems but 
with asymmetric periodic potentials, the ratchet effect has been investigated recently 
by solving the Langevin equation numerically [9]. In the present work we have used 
stochastic form of the 4*^* order Runge-Kutta method. Numerically, one can calculate 
position dispersion at various times and also the distribution of instantaneous velocities 
for the entire period of motion. The drift velocity can be obtained either by taking 
the mean position at a large time and dividing by the time: < v >= limt_*oo ^'^^/^'^ 
or < t; >= vp{v)dv, where p{v) is the velocity distribution. The < ... > in 
< V >= limt^oo "^^^/^^ is the average evaluated over several realizations of the stochas- 
tic run whereas in the later case the (instaneous) velocity distribution is obtained for 
several realizations of the entire stochastic run. In the following we present results 
of our calculation first for the adiabatic case and then for the time-varying (finite 
frequency) field case. 
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FIG. 2: Shows the ratchet current < Av > as a function of F using MCFM (f joined by 
dashed lines) and the Langevin simulation method (* joined by dashed lines) for 7 = 0.035, 
T = 0.4, (p = 0.35. The dotted line for zero current is shown just to guide the eye. 
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FIG. 3: The variation of < > as a function of F is shown for 7 = 0.04 (o), 0.1 (+), 0.2 
(x), 0.5 (open square), 0.75 (*) and 1.0 (solid squares) 




FIG. 5: The ratchet current < Av > when the system was driven by a square- wave field of 
ampUtude Fq are shown (+) with errorbars for 70 = 0.035 and T = 0.4. The curve (with 
open box) is for the adiabatic case for the same 70 and T 
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FIG. 6: The position dispersions < (Ax)^ > are plotted as a function of time t for Fq = 0.08 
(solid line), 0.3 (dashed line) and 0.6 (small dashed line) for 7 = 0.035 and T = 0.4. The 
straight lines are drawn to indicate the average slope equal to 1 and intercepts 1.8, 3.2 and 
2.27 respectively for the above Fq values of the curves at large t limit. 

III. NUMERICAL RESULTS 

A. Adiabatic Case 

The matrix continued fraction method is used to calculate the drift velocity, (and 
hence mobilities), the average potential and kinetic energies, etc, as a function of 70, -fo 
and T. However, in this paper we present only the drift velocities to illustrate the 
ratchet effect in inhomogeneous (non-uniform friction) inertial systems. As will be de- 

n 

tailed in Ref [lOy the expansion coefficients C„(a;) are expressed in terms of their Fourier 
components which form column matrices Cn for each n. By a clever transformation 
Cn and Cn-i are related. Ci being constant it can be found through the normaliza- 
tion condition, 1 = J^^ Co{x)dx and by evaluating the relational matrix (between Co 
and Ci) which depends on the Fourier componets of V{x) and 7(x) and F. A typical 
mean velocity evaluation may require as large as 1000 and the number of Fourier 
components typically as large as 30, if at all a stable solution is to be found. The drift 
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FIG. 7: The position dispersions < (Ax)^ > as a function of time t are shown for the square 
wave drive amplitude Fq = 0.2 but with periods 500 (long dashed line), 2000 (small dashed 
line), 4000 (dotted line), and 8000 (dashed-dotted line) for 7 = 0.035, T = 0.4, and = 0.35. 
The dispersions when a constant force Fq = 0.2 is applied is shown by the solid line. This 
curve is analogous to the homogeneous case of [7]. 

velocity, obtained by taking the sum of mean velocities corresponding to the applied 
force F — \Fo\ and F — —\Fo\, as a function of F is shown in Fig. 2 for 70 = 0.035. 
This case of calculating the drift velocities or the ratchet currents, corresponding to 
the effect of application of a zero averaged drive field of frequency a; — 0, is referred to 
as the adiabatic case. In the same figure we have also plotted drift velocities obtained 
by numerically solving the Langevin equation. Here we have presented the sum Av 
of v for F = ±|Fo|, where v — ]im.t-,oo ^ without the ensemble averaging. The en- 
semble averaging, however, does not qualitatively change the nature of the graph. The 
comparison of the two graphs shows that the Langevin simulation yields qualitatively 
similar result to what is obtained using MCFM. However, the peak ratchet currents 
occur at slightly different values of F, both being at F <^ 1. Importantly, the ratchet 
current, in this adiabatic limit, results only in a very small range of F. In the following, 
for the adiabatic case, we present results obtained from MCFM calculations only. 



From these calculations we find that the ratchet current do not change its qualita- 
tive behaviour for different 70. However, the narrow peaks occur at values of F that 
depend on the values of 70 and T. As 70 is increased the peak position gradually shifts 
towards larger F and also the peak width increases. The largest peak for T = 0.4, 
occurs for 70 = 0.4. Also, at T = 0.4 the ratchet current is in the negative direction 
(Fig. 3). However, there is a distinct possibility of current reversals as a function of 
temperature T as is seen in Fig. 4. The figure also shows clearly that the net current 
could also peak as a function of the temperature T. It is worth reminding that the 
current < Av > here (in the adiabatic case) is just the sum of the average currents 
corresponding to F = ±|-F|. The average currents corresponding to ±|-F| need not, 
however, be nonmonotonous individually. From Fig.3, there is a very clear indication 
that < Af (7) > will also be non- monotonous. 

Since there is no way we can calculate analytically the average currents for time 
varying fields we present (numerical) Langevin simulation results for the non-adiabatic 
case. 



B. Non-adiabatic case 



In the finite (albeit small) frequency drive case, we consider a square wave force field 
where F{t) = Fq for time interval Tq = 1000 and in the next {Tn = 1000) intervalQ 
the force is taken as —Fq and repeated the procedure for a large number of periods 
(typically equal to 5,000 complete cycles) and the mean velocity is calculated. An 
average of this is taken over several realizations (typically 20), thus evaluating the 
ensemble averaged mean velocity or the net (ratchet) current. The obtained steady 
state current as a function of Fq is plotted in Fig.5 for 70 = 0.035. For comparison 
the current obtained in the adiabatic case is plotted together. The range of F for 
which ratchet current can be obtained has broadened considerably in comparison to 
the adiabatic case. The value of F at which the peak current is obtained has also 



shifted to a higher value. In Fig. 6, one can see that around t ~ 10^, the position 
dispersion < (Ax)^ >=< x'^(t) > — < x{t) t hinting at the approach to the 

diffusive steady state situation. 

Fig.6 with < {Ax{t)y >~ t provides a method to calculate the diffusion constant D 
in the large time limit. It is interesting to observe that in this limit < (A,t)^ > shows 
a peaking behaviour as a function of F. For 7 = 0.035, the peak ratchet current occurs 
close to where < (Aa;)^ > is maximum. It would, therefore, appear that ratchet current 
in this case is dispersive and not coherent. However, when the relative dispersion is 
plotted as a function of F, ^^^^^2^ shows small values < 1 in the range of F where 
ratchet current is obtained. Thus, the obtained ratchet current is, indeed, a result of 
coherent motion of particles. Equivalently, the related Peclet number Pe which is a 
measure of coherence of motion, turns out to be much larger than 2 indicating the 
particle motion giving rise to steady state ratchet current to be coherent. The motion 
at the intermediate (transient) time scales, however, shows a very different nature. 

Q 

As mentioned earlier, Lindenberg, et. al.'Tl] have shown that the motion of an 
ensemble of non-interacting particles in a periodic potential together with a uniformly 
applied constant force showed dispersionless behaviour at intermediate time scales. 
They have set the time interval [tmin < t < tmax] for a given applied force Fq during 
which the system shows such an interesting behaviour. This behaviour is seen only 
for a certain range of Fq values. A typical graph of < (Ax(t))^ > on log-log scale is 
presented numerically in Fig. 7. In this figure are also shown the effect of frequency of 
the dichotomously time varying (square wave) applied field of same amplitude Fq = 0.2. 
If the half period of the applied field over which the field remains constant before its sign 
is changed, Tq < tmin, then after Tq the dispersion shows a dip before again picking up. 
This is repeated in every Tq. If, however, Tq > tmin a.t t = tmin the slope of < Ax^ > 
changes abruptly. At t = Tq, < (Ax)^ > begins to increase sharply before it appears 
to level itself up till it reaches t = 2Tq. The behaviour shown between Tq and 2Tq is 
seen to repeat every Tq for a long time. The frequency dependence of applied force on 



the dispersive behaviour of particle motion shows interesting possibihties of dispersion 
reduction as well as brief dispersionless motion. 



IV. DISCUSSION AND CONCLUSION 

The present work gives numerical evidence of the possibility [4] of obtaining ratchet 
current in an inertial inhomogeneous symmetric periodic potential system subjected to 
a symmetrically driven periodic external field. The nonuniform temperature case has 
been investigated earlier by Blanter, et al. |l2|. We have calculated the net current from 
the average distance travelled after a long time. However, almost the same current is 
obtained from the (instantaneous) particle- velocity distribution over the entire sojourn 
time of the particles. But we find that the efficiency with which useful work can 
be derived from such a system is quite low and falls in the subpercentage range. This 
could be because the system inhomogeneity considered here is a feeble cause to generate 
net current without the asymmetric support of external forcings. Preliminary work, 
however, shows that the magnitude of ratchet current can be increased dramatically if, 
instead of time symmetric drive, time asymmetric drive with zero time averaged force 
per period is used. Similar effect for overdamped case has been reported earlier jisj]. 

The velocity distribution in this case shows an overwhelming asymmetry in favour 
of particle motion in the direction of current in contrast to the symmetric drive case 
where it shows almost symmetric velocity distribution. The width of the velocity 
distribution 9| increases with the amplitude of the periodic driving. In conclusion, 
even though system inhomogeneity is a weak ingredient to obtain ratchet current it is 
a distinctive way of obtaining such a current. 
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